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A brief review is given of recent developments in the study of surface critical phe- 
nomena from the viewpoint of exactly solved lattice models. These developments 
include exact results for the polymer adsorption transition and the surface critical 
exponents of the eight-vertex model. 



1 Introduction 

There can be little doubt that our understanding of phase transitions and 
critical phenomena has been greatly enhanced by the study of exactly solved 
lattice models in statistical mechanics. 1 The early pioneering work of Bethe, 
Onsager, McGuire, Yang, Lieb and Baxter has left us with the legacy of the 
Yang-Baxter equation, which governs integrability in the bulk. Many models 
have been investigated to yield bulk critical behaviour. E.g., Baxter 1 ' 2 derived 
the bulk free energy of the 8-vertex model, and thus the bulk specific heat 



where t = T — T c is a measure of the deviation from criticality and ctf, = 2 — ir//i 
is the continuously varying bulk specific heat exponent, with ctf, = O(log) in 
the Ising limit. 

In general there is an i?-matrix R(u) satisfying 



whose elements define the vertex or Boltzmann weights of exactly solvable 
lattice models in 2-d statistical mechanics. On the other hand, Sklyanin has 
given a formulation of integrability at a surface, governed by both the Yang- 
Baxter equation and the reflection equation, which reads 3 ' 4 

R 12 (u - v)K^ (u)R 21 (u + v)K^ (y) = (v)R 12 (u + v)K^ (u)R 21 (u - v) (3) 

where K~ (u) is a reflection or 7f-matrix for the right boundary. A similar 
equation holds for K + (u) at the left boundary, which is automatically satisfied 
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as t -> 



(1) 



Ri 2 (u)R 13 (u + v)R 23 (v) = R 23 (v)R 13 (u + v)R l2 {u) 



(2) 
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if K + (u) = K~ (— u— 7j)'M, where r\ is a crossing parameter and M is a crossing 
matrix (both being specific to a given i?-matrix). 4 ' 5 

Integrable quantum spin chains with open boundaries, such as the XXZ 
chain, 6 are recovered naturally in the language of 7f -matrices. 4 ' 7 Most of the 
attention on spin chains has been on the quantum group invariant cases. 7 ' 8 
A'-matrices have been found for a number of models, although up until quite 
recently very little in the way of surface critical phenomena has been derived 
from them. 

The most natural open geometry for vertex models is that depicted in 
Fig. 1. The integrable boundary weights follow directly from the Tf-matrix 
elements. Specifically, the left and right boundary weights are given by 9 

Z»(«) = 5>£(2««(«), r b a (u) = K~ b (u). (4) 

ab 

I give a brief account in Sec. 2 of results obtained with Ming Yung for the 
dilute 0(n) loop model with open boundaries and related self- avoiding polymer 
problems. 8-12 I then discuss more recent work 14 with Yu-Kui Zhou in Sec. 3 
on surface critical phenomena in the 8-vertex model." Further developments 
are mentioned in Sec. 4. 




Figure 1: The open square and honeycomb lattices. 



2 The Dilute O(n) Loop Model and Related Polymer Problems 

The dilute O(n) loop model on the square lattice 15 turns out to be solvable 
for two sets of boundary weights. 9 The partition function is defined as 

^oo P = ^pr---pr 3 i3 " p (5) 

g 

a Our results for the 8-vertex model were announced at this Meeting. 
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where the sum is over all configurations Q of non- intersecting closed loops which 
cover some (or none) of the edges of the square lattice in Fig. 1. The possible 
vertex configurations are shown in Fig. 2, with a vertex of type i carrying a 
Boltzmann weight p t . In configuration Q, m l is the number of occurrences of 
the vertex of type i, while P is the total number of closed loops of fugacity n. 




Figure 2: The allowed loop model vertices and corresponding weights. 

The integrable bulk loop weights are 15 ' 16 

pi = sin(3A — u) sin(u) + sin(2A) sin(3A) 

P2 = p3 = ei sin(3A — u) sin(2A) 

Pi = P5 = £2 sin(u) sin(2A) 

P6 = P7 = sin(3A — u) sin(u) (6) 

ps = sin(3A — u) sin(2A — u) 

Pq = — sin(A — u) sin(u). 

where e\ = e| = 1 and n = — 2cos(4A). The two sets of integrable boundary 
loop weights are 9 

Pio = Pi3 = sin[i(3A - u)] 

(A) i (7) 
Pn = P12 = ei sm[ 5 (3A + u)] 

Pw = Pi3 = cos[|(3A - u)] 

(B) 1 (8) 
Pn = P12 = ei cos[2 (3A + u)]. 

These weights follow from the A'-matrices of the underlying Izergin-Korepin 
vertex model. 8 The eigenvalues of the transfer matrix have been determined 
for both cases by means of the Bethe Ansatz. 9 



Exact results have been obtained for various polymer problems on the 
related honeycomb lattice, which follows in the honeycomb limit 15 u = A 
(where the bulk weight pg =0). In this way loop configurations on the square 
lattice become loop configurations on the honeycomb lattice depicted in Fig. 
1, where the partition function reduces to 10 



J \oop 



J2* L y L °n p . (9) 



At n = the contribution from all closed loops vanishes and the partition sum 
reduces to a generating function for self-avoiding walks where x is the fugacity 
of a step in the bulk and y is the fugacity of a step along the surface, in this 
case either edge of the strip. Here L is the length of the walk in the bulk and 
L s is the length of a walk along the surface. 

The values of x and y are fixed by integrability. The well known Nienhuis 
result 17 



l/x* = y2±V2 Zr n (10) 

follows from the bulk weights. On the other hand, at the boundary y* = x* 
for case (A), i.e., the critical surface coupling is equal to the bulk coupling. In 
the language of surface critical phemonena, 18 this corresponds to an integrable 
point on the ordinary transition line. For case (B) 

1/y* = \l±V2^n~. (11) 

This point corresponds to the 0(n) model at the special transition. 

The polymer adsorption transition corresponds to the n = limit of the 
O(n) model at the special surface transition. 18 At n = we have 



1/x* = V / 2~+\/2~ and 1/y* = \f^2. (12) 
The exact critical adsorption temperature T a is thus given by 10 ' 11 

exp (j^j = y*/x* = yJl + y/2 = 1.553 . . . (13) 

where e is the contact energy at the surface. Physically, there is a desorbed 
phase for T > T a and an adsorbed phase for T < T a . The number of configu- 
rations of self-avoiding walks with one end attached to the boundary has the 
asymptotic form 19 

/ fi L L^ T > T a 

^L-rr-i T = T a (14) 
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where for the honeycomb lattice, /j = 1/x* . The surface critical exponents 7° 
and 7^ p can be obtained from the finite-size corrections to the transfer matrix 
eigenvalues, with result 10 ' 13 ' 20 

7l ° = M and 7l sp = f. (15) 

The value for 7° is the conformal invariance result of Cardy, 21 while the value 
for 7j P had originally been conjectured by Guim and Burkhardt. 19 The exact 
value of the crossover exponent also follows as <f> = \ . 
A further configurational exponent, 

7i = I, (16) 

is obtained by considering mixed ordinary and special boundary conditions 
arising from different 7f -matrices on either side of the strip. 12 ' 13 This case can 
also be solved exactly by means of the Bethe Ansatz and the exponent 71 
is seen to control the number of self-avoiding configurations beginning from 
near some origin at the surface with a non-adsorbing boundary to the left 
and an adsorbing boundary to the right (see Fig. 3). The same adsorption 
temperature T a holds for the mixed boundary case. These results have been 
numerically confirmed in the half-plane using series expansion techniques. 22 

Although we have obtained various exact exponents by explicit calculations 
on the honeycomb lattice, they are expected to be valid for other planar lattices 
on universality grounds. 



Figure 3: A self-avoiding walk from an origin on the surface of the half-plane. 



3 The 8-vertex Model 

A'-matrices have also been found for the 8-vertex model. 23 As for bulk weights 24 ' 
the boundary weights can be interpreted in terms of Ising interactions (see Fig. 
4) . The parametrisation of the vertical and horizontal Ising couplings K and 
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Figure 4: The bulk and surface vertex and Ising spin configurations. 



L and the four-spin interaction M is given in terms of elliptic theta functions 



, 24,25,1 



ac 
bd 



l(«)J 



„4L 



ad 



l(A-u) 



,4M 



l(A) 



Recalling (4) , the Ising coupling A' s in the surface layer follows in terms of the 
known A'-matrix as 



,4/4' 



».2„1 



K u (u/2)K 22 (u/2) 



(17) 



A'r 2 (w/2)A' 2 - 1 ( U /2) 

In this way we can realize an anisotropic Ising model with four-spin in- 
teractions with open boundary conditions (see Fig. 5). Although we have the 
integrable boundary weights, so far there is no exact solution for this model. 
However, an inversion relation for the free energy has been obtained by Yu-Kui 
Zhou. 26 We have used the inversion relation method 1 to solve this relation for 
the surface free energy and thus the surface specific heats 18 C s and C\ . We 
find that as t — > 0, 



t 



t 



with 



with 



2- — 

2,j/ 



(18) 
(19) 



The known Ising results, 18 a s = l(log) and a\ = — l(log), are recovered at 
(j, = it/2. Our results are in agreement with the scaling relations a s = <Xb + v 
and a.\ = ah — l, 18 where v = ir/2/i is the correlation length exponent. 27 
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Figure 5: The 8- vertex lattice (dotted lines) and the Ising lattice (broken and solid lines). 

4 More Recent Developments 

Sklyanin's approach can be reformulated in the language of face models in 
which the A'-matrices are interpreted as boundary face weights of a double row 
transfer matrix. This has recently been done by a number of authors. 26 ' 28-31 
Such A-matrices were presented at this Meeting for the Andrews-Baxter- 
Forrester model in the talks by Koo and Pearce. 30 ' 32 There have been some 
rapid developments in this area since this Meeting. First, we repeated our 
8-vertex model calculations for the Andrews-Baxter-Forrester model to derive 
the excess specific heat exponent a s = (7 — L)/4, with the Ising value again 
recovered for L = 3. 33 We also derived the exponents a s and a\ for the CSOS 
model 34 and the exponents 5 S and a s for the dilute Al models. This latter 
work yielded the value 5 S = — ¥^ f° r an Ising model in a magnetic field. 35 Fur- 
ther A'-matrices have been found for general face and vertex models associated 
with the algebras A™ , , CP , and A^? 6 
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